Theory of Thermal Motion in Electromagnetically Induced Transparency: 
Diffusion, Doppler, Dicke and Ramsey 
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We present a theoretical model for electromagnetically induced transparency (EIT) in vapor, 
that incorporates atomic motion and velocity-changing collisions into the dynamics of the density- 
matrix distribution. Within a unified formalism we demonstrate various motional effects, known 
for EIT in vapor: Doppler-broadening of the absorption spectrum; Dicke-narrowing and time-of- 
flight broadening of the transmission window for a finite-sized probe; Diffusion of atomic coherence 
during storage of light and diffusion of the light-matter excitation during slow-light propagation; 
and Ramsey-narrowing of the spectrum for a probe and pump beams of finite-size. 



I. INTRODUCTION 

The Doppler effect, discovered in the mid-19th cen- 
tury, causes a broadening of spectral lines in thermal 
media which is linearly proportional to the radiation 
wave- vector In 1953 R. H. Dicke predicted that 

the Doppler-broadened spectrum can be narrowed due 
to frequent velocity-changing collisions 0, 0, 0], by a 
factor proportional to the ratio between the collisions 
mean free-path and the radiation wavelength. This phe- 
nomenon, known as Dicke narrowing .was observed for 
microwave and optical transitions [5j, [fj . When the mo- 
tion of the atoms is diffusive, the resulting width is pro- 
portional to the diffusion coefficient and quadratic in the 
radiation wave- vector 0, 0]- Therefore it is sometimes 
referred to as diffusion-narrowing (of the Doppler pro- 
file) or diffusion-broadening (of the spectrum of a sta- 
tionary atom). For a finite-size beam, as illustrated in 
Fig. [U both the Doppler and the Dicke widths can be 
explained as a time-of-flight (TOF) broadening. A com- 
prehensive literature survey and a theoretical treatment 
of the Doppler-Dicke problem is presented by May @. 

Electromagnetically Induced Transparency (EIT) oc- 
curs when two light fields, a probe and a pump, couple 
two sublevels of an atomic ground state manifold to a 
common upper level. When the Raman resonance con- 
ditions are satisfied, a quantum coherence evolves within 
the ground state in the process of coherent population 
trapping (CPT) [9j , inducing a unique transparency win- 
dow in the absorption spectrum, accompanied by dis- 
persion. The spectral width of these features depends 
primarily on the ground-state decoherence rates, ranging 
from the order of several Hz in cryogenically-cooled crys- 
tals [l(| up to the order of several KHz or MHz in ther- 
mal vapor pj| . A wide variety of phenomena has been 
demonstrated utilizing the ultra-narrow EIT resonances, 
e.g. slow light [12], stored light [131 ] . and non-linear op- 
tics at low light levels 14]. Having the advantage of 
being relatively simple to implement, thermal vapor EIT 
has been used for various applications, such as frequency 
standards [H, EH, magnetometers [TtJ and optical de- 
lay lines. In this work we present a comprehensive ana- 
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FIG. 1: Illustration of time-of-flight (TOF) broadening in the 
Doppler (left) and Dicke (right) limits. Assume a beam of 
width Ax in the transverse plane and width Ak± ~ 1/Ax 
in fc— space. Atoms with a transverse velocity v t h (left) cross 
the beam in time Aa; /v t h and cause a TOF broadening of the 
order of v t hAk±, which is equal to the well-known Doppler- 
width. Atoms that undergo diffusion (right), traverse the 
beam in average time of Ax 2 /D, where D — Vt^A is the diffu- 
sion coefficient and A is the mean-free path between collisions. 
This results in a TOF broadening of the order of vth^Ak\, 
which is the well-known Dicke-width. 



lytic model for the effect of atomic motion on vapor EIT, 
describing broadening and narrowing mechanisms, and 
the phenomena of diffusion during slowing and storage 
of light, within a unified formalism. 

As a two-photon process, EIT in room-temperature 
vapor is susceptible to residual Doppler broadening, pro- 
portional to Sq— | q x — q2 1 , where qi and q2 are the wave- 
vectors of the probe and the pump, respectively. However 
in the presence of buffer gas, when diffusion dominates 
the atomic motion, a Dicke-like narrowing of the Doppler 
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spectrum may occur and the narrowing factor is propor- 
tional to the ratio between the mean- free path and Sq _1 
[HI]. In an hyperfme EIT, Sq^ 1 is of the order of cen- 
timeters and the narrowing factor is of the order 10~ 4 , 
completely diminishing the residual Doppler effect and 
allowing the high accuracy of EIT-based frequency stan- 
dards (la . [T^ | . By introducing a small angular deviation 
between the probe and the pump, it is possible to quan- 
titatively measure the residual-Doppler width and the 
Dicke width of Zeeman EIT lines [2(1 [HI . 

For probe and pump beams of finite size, the EIT spec- 
trum is subjected also to TOF broadening [22J. Here 
again, the broadening reduces to the Dicke-type in the 
limit of diffusional motion. Nevertheless, the spectrum 
is also affected by a more complicated mechanism, re- 
cently denoted as Ramsey narrowing, which is attributed 
to the random-walk of atoms in and out of the beams 
[23L l24j ] . The process in which an atom carrying ground- 
state coherence leaves the beam, evolves "in the dark", 
and then re-enter the beam, is analogous to the Ram- 
sey method of separated oscillating fields 25]. This pro- 
cess is more relevant for two-photon phenomena, such as 
velocity-selective optical pumping, magneto-optical spec- 
troscopy and CPT/EIT, because of the small frequencies 
(c ■ Sq) and long coherence times, that are comparable 
with the typical time the atoms spend inside and out- 
side the beam. Ramsey-like features in such systems 
were observed in the velocity, time and space domains 
[26L I27L [H, Ramsey narrowing occurs due to the 

random distribution of the durations " in the dark" , re- 
sulting in a superposition of the spectral Ramsey fringes, 
that wipes out the fringes and leaves a single narrow fea- 
ture in the center. A theoretical technique of integrating 
over all possible Ramsey sequences, taking into account 
calculated distributions of the durations in and out of the 
beams, agrees well with experimental results [241 ] . 

While Ramsey narrowing is an evidence for the dif- 
fusion of EIT coherence, a more explicit demonstration 
was given in recent storage-of-light experiments, in which 
an elaborated spatial profile of the probe field was stored 
and retrieved in vapor [3(1 [HI]. The restored profiles were 
exactly predicted by assuming that during storage, when 
no fields are present, the ground-state coherence under- 
goes regular diffusion. Since the coherence is complex, 
the phase pattern also diffused and interference effects 
were observed. Neither a more basic theory that derives 
the diffusion equation for the atomic coherence in the 
absence of fields nor a prediction for the diffusion-like 
behavior of slow-light were presented. 

The model presented in this paper describes the ef- 
fect of velocity-changing collisions on the position and 
velocity dependence of the atomic density matrix. As- 
suming a Boltzmann-like relaxation in velocity space, we 
write in section[TT]the dynamics equations for the internal 
and the external atomic motion, under the weak-probe 
approximation. We then derive the equations of motion 
for the slowly- varying envelopes of the atomic coherences 
and the probe's field, while the pump's envelope is as- 



sumed stationary. In section HTT1 we formally solve these 
equations and derive the probe's susceptibility, for the 
general case and for a plane-wave pump. By general we 
mean that neither the "Doppler" nor the "Dicke" lim- 
its are taken [8( for the one-photon or the two-photon 
transitions. This model extends our previous study on 
Dicke-narrowing [l8j in the following aspects: it avoids 
the use of the low-contrast approximation, incorporating 
the power-broadening effect; it allows for a non- planar 
(finite-size) probe beam, introducing a wave-vector filter; 
it is time-dependent and thus allows for the propagation 
of a probe package. We note also that the assumption 
of Gaussian distribution of atomic trajectories is avoided 
here (see Eq. (9) in Ref. [H|). In section ITVl we consider 
the realistic regime, where the velocity relaxation-rate is 
large enough to cause Dicke-narrowing of the two-photon 
transition. In this regime we derive a diffusion equation 
for the density-matrix distribution during storage-of-light 
and afterwards calculate the dynamics in the presence of 
the fields. For the case of plane-wave pump we recover 
the Dicke-Diffusion absorption spectrum and arrive at a 
diffusion-diffraction equation for the slowly propagating 
envelope of the probe. For the case of finite-size pump 
and probe, analyzed in section [Vj we solve the diffusion 
equations and retrieve the Ramsey-narrowed absorption 
spectrum. Note that the latter is done in a steady- 
state approach, rather then by averaging over all possible 
atomic trajectories. 

II. EQUATIONS OF MOTION 

A. Atom-Field Interaction 

We consider three-level atoms in a A— configuration, 
with an upper state |3) and two lower states, |2) and 
|1) , as depicted in Fig. [3J The atoms interact with two 
external, classical electromagnetic fields, propagating in 
time t and space r, 

E,(r,i) =sfte{— e s Q s (r,i)j far a =1,2, (1) 

with 

n s (r,f)=O s (r,t)e-^V^' r , (2) 

where e s and q s are the polarization vector and wave 
vector of the probe (s = 1) and the pump (s = 2); 
lu s = c|q s |; /i3 S = (3|d-e s |s) is the 3 — ► s transition 
dipole moment matrix element; and fl s (r, t) is the slowly 
varying envelope in time of the Rabi frequency, satisfying 
|f O s (r,i)| <w s |Q s (r,t)|. 

First we consider a single atomic system, the i— th 
atom, one of many identical particles. The center of 
mass position rj (t) enters the internal dynamics of the 
atom because the external fields must be evaluated there. 
We denote the ss' reduced density-matrix element of the 
i— th atom as p l ss , (f) . For a weak probe, namely when 
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FIG. 2: (a) A probe beam and a pump beam, with a finite 
envelope in space, propagate through the vapor cell. The 
z— axis is chosen perpendicular to the probe's direction and 
the x and y axes form the transverse plane, (b) Atomic levels 
diagram. Sli and 0.2 are the Rabi frequencies of the probe 
and the pump, respectively. Td and Ti2 are the decoherence 
rates of the optical and the ground-state transitions. 



where g = 2iruj\ |/i3i| 2 /H is a coupling constant and V 2 
is the three-dimensional Laplacian operator. To obtain 
Eq.© we neglected the second order temporal deriva- 
tives of the envelopes and chose, without loss of general- 
ity, qi = zqi, where z is a unit vector in the z— direction 
(see Fig. [5]a). The second order spatial derivatives are 
retained to allow the description of finite-size beams. 



B. Boltzmann-like Equations 

We consider a dilute gas of active atoms in the pres- 
ence of a noble buffer-gas that causes frequent velocity- 
changing coherence-preserving collisions. We introduce a 
generalized density-matrix distribution function in space 
and velocity 

Pas' = Pes' (r, V,t) = X; A.' (*) S(T-Ti (t)) 5 (v - V, (*)) , 

i 

(?) 

where the time dependence of p\ s , (t) is due to the in- 
ternal dynamics only, considered earlier in the single- 
particle density-matrix. We express 



the probe field is weak enough to cause p\ t , the 

equations of motion of /9 31 and p l 21 are (cf. jHI) 



T d p z 31 - iuJ 3 iPli + (Ti,t) p\ x 



Q . 

-Q- t P21 = ~ r 2lP21 ~ i^2lP\l + l ^2 ( r i,t) Plf 



(3) 



Here W31 and are the frequency and the decoherence 
rate of the 3 — > 1 transition; LO21 and T21 are the fre- 
quency and the decoherence rate of the 2 — > 1 transition; 
p l 31 is linear in the probe field; and the equilibrium den- 
sity matrix in the absence of the probe (Oi = 0) is 



d 

dl Ps 



E 



|/4'MU- -.i-,(/!M(v-v,(M! 



Pss> 



^^(r-r,(t))^. 



dl ps 



E 



dt 



dl Ps 



^(v-v i( t)) 



col . 



-pl s/ )S(r-v t (t))5(v-v t (t)), (8) 



s = s' = 1 
otherwise. 



(4) 



We consider the case of non-depleted pump, i.e. we 
assume that the pump's envelope is constant in time and 
is given by CI2 (r, t) = (I2 (r) . The wave equation of the 
probe field is 

, 1 d 2 



dt 2 



4tt d 2 

E i( r ^) = 5^2 p 3l(r,t), (5) 



where P 31 (r, t) = 5Re |p 31 (r, t) e -iwiV qi ""} is the con _ 

tribution of the 3 — > 1 transition to the expectation value 
of the polarization density, and P31 (r, t) is the slowly 
varying polarization. With Eq.fTJ) we have 



d 







dt dz 



2qi 



fJ 



^l+ c 7T- l TT~ V "1 (r, t) = i^-ei ■ P31 (r, t) 



M31 



(6) 



where ^ 



is the effect of collisions. Notice that 



the density-matrix distribution function is classical, as 
far as its external motion is concerned, and quantum 
mechanical in its internal atomic motion. The function 
p SS ' (r, v, t)d 3 rd 3 v is interpreted as the probability to find 
an atom of density-matrix element p ss i at the time t, near 
position r, within a volume element d 3 r, with velocity v, 
within a velocity volume element d 3 v. 

With this interpretation we understand the effect of 
collisions as relaxation to thermal equilibrium of the cen- 
ter of mass, and we shall take it into account using the 
Boltzmann relaxation method. Introducing the density 



R ss > (r, t) = J d 3 vp ss , (r , v, t) , 



(9) 



which is interpreted as the number of atoms with p ss > 
per unit volume near r in space, we assume that the 
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distribution p ss '(r, v, t) relaxes to equilibrium as 



col . 



p ss> (r, v, t) = -7 [p ss > (r, v, t) - R ss > (r, t)F (v)] . 

(10) 

Here 7 is the relaxation rate in time, proportional to 
the collisions rate, and F (v) is the thermal equilibrium 
Boltzmann distribution function in velocity space, 



F(v) 



1 



(2-4) 3/2 



-« 2 /(2< h ) 



. (11) 



where kg is the Boltzmann constant, T is the tempera- 
ture, m is the atomic mass and Vth is the thermal velocity. 
Here the internal motion, and the external one, are com- 
pletely separated. The collisions with the buffer gas af- 
fect only the external degrees of freedom and are assumed 
much frequent than the collisions with the active atoms. 
If the buffer gas particles affect the internal motion, e.g. 
by pressure-broadening Q , it is taken into account in the 
atomic decay rates. The collision term, Eq. lfTTJ)) . neglects 
partial wave scattering for I > and also neglects the 
energy dependence of the scattering cross-section. From 
Eqs. (|3l l4l8[) and the Boltzmann collision term of Eq. (fl"0]) . 
we observe that 

£ p^ cq) S (r - r, (t)) 6 (v - v< (t)) = n F (v) , (12) 

i 

where no is the gas particle density, and the equations of 
motion are: 

d d\ 

— +r d + zw 3 i +v-— ) psi - ifl 2 {r)p2i (13) 
at or J 

(r, t)n Q F (v) + 7 [p 31 - R 31 (r, t)F (v)] = 0, 

d d\ 
-Q- t + r 2 i + 1U21 + v-— J p 2 i - i£l* 2 {r)pzi 

+7 [p2i -i?2i(r,t)F(v)] =0, 

with P21 = p2i(r,v,t) and p 31 = p 31 (r,v,t). We note 
that the above semiclas steal analysis of the dynamics 
has a quantum-mechanical equivalent, which will be of 
significance in a regime where the atomic recoil velocity 
is comparable with the thermal velocity. 



C. Envelope Equations 

Since the probe field propagates through the cell with 
a slowly varying envelope, p 3 i(r, v, t) and p 2 i(r, v, t) can 
be expressed as 



p 31 (r,v,i)=p 31 (r,v,t) e -^ lt e lqi - r , 
pai(r,v,t) = p 21 (r,v,i) e -^-^V^-^- r , 



(14) 



where p 3 i and p 2 \ are slowly varying in space and 
time. Similarly we introduce the slowly varying densi- 
ties, i?2i(r,t) and i? 3 i(r, t), and express the expectation 



value of the polarization density, P 3 i (r, f), in terms of 
the density -ff 3 i (r, t) as 



ei -P 3 i (v,t) = p,* 31 R 31 (r,t). 



(15) 



With the one-photon detuning, Ai = u>i — w 3 i, and the 
two-photon Raman detuning, A = u>±— u; 2 — w 2 i, we define 



$1 = A 1 -q 1 -v + i(r d + 7 ) ) (16a) 
e 2 = A-( qi -q 2 )-v + i(r2i+7). (16b) 



and write Eqs. fUS)) and Eq.© as 



di+v-fr-iti) p3i(r,v,t)- 7j R 31 (r,t)F(v) (17a) 

-jf2 2 (r)p 21 (r,v,i) - ffli(r,t)n F(v) = 0, 
di+v-fr-ib) p2i(r,v,i)- 7j R 2 iM)F(v) (17b) 



-in5(r)p 3 i(r,v,t) = 0, 



and 



» +c |- i ^-V a )fi 1 (r f *) = ipft 1 (r f *). (18) 

Eas. (fT7|) and (JT5J) compose the full set of equations of 
motion for the slowly varying envelopes. 

Finally, in sections IIIII and IIVI we study the case of 
a stationary, plane- wave pump, fi 2 (r) = ^2- For this 
case it is convenient introduce the Fourier transforma- 
tion and replace the slowly varying time dependent and 
r— dependent functions by 



/M) 



and write Eqs. (|T7|) as 

(w - k • v + & ) psi (k, v, u) - i~/R 3 i (k, uj) F (v) (20a) 

+Sl 2 P2i(k,v,w) +fii(k,w)n F(v) = 0, 
(w - k • v + £2) P2i(k, v, lu) - i 7j R 21 (k, lu)F (v) (20b) 

+ft* / 5 31 (k,v,^) =0, 

and Ea.([T8]) as 



ikz _ ^ +i ) Oj (k,w) = i^flai(k,w). (21) 
c zqi / c 



III. GENERAL SOLUTION 

We consider the EIT medium in a cell, with a probe 
beam of finite width in the transverse plane (x,y), prop- 
agating along the z— axis in the direction of qi = zqi. 
The variation of f^i in space is much slower than 2n/qi 
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and the variation in time is much slower than 2~k/u>i_. 
The pump is a plane- wave, 2 , propagating with a wave- 
vector q 2 . Here we present the general solution for the 
probe field propagation inside the cell, i.e. away from its 
boundaries, for any relaxation rate 7, without taking the 
"Doppler" or the "Dicke" limits. 

We start from Eqs. (|2"0]) and formally solve for p^i and 

P2l- 



Pai = 
and 
hi 



Fiy) 

F(v) 



ilRsi ~ fii(k, u/)n 2 

17^21 w — kv + £ 2 



, (22) 



oj - k • v + £1 ijR 3 i — fii(k, oj)n 



ijR 2 



(23) 



where | • • • | stands for matrix-determinant and 

U = (lj - k • v + a) (w - k- v + &) - M2I 2 . (24) 

Integrating over velocity, we get coupled equations for 
^31 and R 2 i, 





R31 








R21 _ 




f^G 



1 - ijGx i^n 2 G 
i^Cl* 2 G 1 - ijG 2 



where we have introduced the integrals 
G(k,w) = / d 3 «F(v)l, 
Gi(k,w) = /" d 3 vF(v) 
G 2 (k,w) = / d 3 vF(v) 



-k- 



a; — k • v + £1 



Qi(k, u>)n , 
(25) 

(26a) 
(26b) 
(26c) 



Solving Eq. ((25|) for -R31 we obtain 



6 n \ on n «o /I - *7G 2 (k,w 
i? 3 i (k,u>) = ifti (k,w) — — — — — - — 
7 V G d (k,w) 



- 1 



where 

G d (k, w) = (1 - ijd) (1 - z 7 G 2 ) + j 2 \n 2 \ 2 G 2 . 



(27) 



(28) 



We can now return to Eq. (|2"Tj) and solve for Cli (k, u>). 
Since R31 is linear in we introduce the linear suscep- 
tibility by, 



i? 3 i (k, w) = X31 (k, w) -fii (k, uj) 
9 



and the complex wave-number, 

■„2 



p(k,w) = h%3i (k,w) 

c 2gi 



to express Eq. (|2Tj) as 

[k z -p(k,w)]fii(k,w) = 0. 



(29) 



(30) 



(31) 



When the changes in the envelopes along the z— direction 
are much smaller compared to the changes in the trans- 
verse plane, we may replace k — > kj_ in the £'s and the 
G's in Eqs. (j2"4"]) and (j2"B"|) , where kj^ is the projection of 
k onto the transverse plane, and write Eq. (|31[) in the 
(z;kj_,w) coordinates: 

d - 

— fii (z;k±,u) = ip(kj_,u))Qi(z;kx,uj), (32) 
oz 

or, the solution 

f2i(z 2 ;kj.,w) =ni(ai;kj.,w)e*( kj -' w ^« a -* I \ (33) 
Finally, the probe's envelope is given by 

J_oo 27T 27T 

x Hi {zi;k ± ,u)e ip( - kx ' u) ( Z2 - Zl \ (34) 

The Doppler-Dicke Transition 

It is first instructive to consider the case of ordinary 
one-photon absorption by taking £^2 = 0, and to assume 
for simplicity a plane- wave stationary probe, namely lj — 
and k = 0. For this case Eqs.(|26" l) - (p9j) give 



X 3i (A 1 )=^JC(A 1 ), 
c 

where K (AO is the one-photon complex spectrum, 

iGi (Ai) 



1-ijGx (AO' 



(35) 



(36) 



and Gi (Ai) takes the form of a Doppler-like profile (a 
Voigt convolution), with 7 being added to the standard 
homogenous- width, r d : 



Gi (AO 



2irv th 



du- 



-«V(2< h ) 



Ai - qiu + i(T d + 7) 



(37) 



The spectrum K(A.±) in the form of Eq. Q36p was previ- 
ously presented by May [811 for one-photon transitions. Its 
extreme limits are 3, 8, 18]: the Doppler limit, trivially 
obtained by setting 7 = 0, and the Dicke limit, found for 
large 7. The Dicke parameter is w t h<?i/7, proportional to 
the ratio between the mean free-path and the radiation 
wavelength. 

For optical transitions in room-temperature vapor, the 
one-photon line is usually in the far Doppler limit, i.e. 
K w iG\. We have defined K such that it is real for 
|Ai| <C \q\u + i (Td + 7)| , i.e. near the one-photon res- 
onance where experiments in EIT are often done, and 
is equal to the on-resonance absorption (in frequency 
units). Specifically, for an atom at rest (7 = and no 
Doppler), K{ A x =0) = Tj 1 . 

A similar Doppler-Dicke transition occurs for the EIT 
transmission window. In order to demonstrate that, 
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FIG. 3: Normalized EIT transmission spectra, numerically 
calculated from Eqs. (f26)l - (f29)l . for three values of the EIT 
Dicke parameter, r\ = Vthk/j (the ratio between the resid- 
ual Doppler width and the velocity relaxation rate), with 
r d = 2500w th fc, r 2 i = 0.025v th k and \Q 2 \ 2 /T d = r 21 /25 
(small power-broadening). When rj is large (solid black line) 
the spectrum is a Voigt curve (a Gaussian-Lorentzian convo- 
lution). When r\ is small (dashed blue line) the spectrum is 
a pure Lorentzian of width Yd- The dot-dashed green line 
demonstrates an intermediate result. 



FIG. 4: Full width at half maximum of the EIT transparency 
window (points), obtained from numerical results similar to 
Fig. [3] as a function of the wave-vector difference, k, for 
various values of the velocity relaxation rate, 7. The dashed 
lines are given by Eq. (|38p . Other parameters are (typical 
for experiments with small power-broadening): v t h — 170 
m/s, T d = 100 MHz, F 2 i = 1 KHz, \n 2 \ 2 /T d = 40 Hz. The 
three lines in Fig. [3] correspond here to \k\ ~ 1.5 mm -1 and 
7 = 16, 160 and 1600 KHz. 



we have chosen a set of typical parameters with small 
power-broadening and calculated several EIT lines from 
Eas. ([25 ]) - ([23 ]) by numerically integrating Eqs. (HU). Three 
line shapes as a function of the normalized Raman de- 
tuning are presented in Figs. [H The full- width at half- 
maximum (FWHM), as a function of |k| and for vari- 
ous values of 7, is presented in Fig. 2J For the cal- 
culations we took a stationary (lo — 0), collinear and 
degenerate (qi = q2) EIT with klz, so that k = |k| 
is the wave-vector difference between the pump and the 
probe (similar results are obtained by replacing k with 
|qi — q 2 |, when |qi — q%\ <C The residual Doppler 
width is expected to be v t hk, and the EIT-Dicke parame- 
ter is r\ = Vthk/^j [HI]. Fig. [3] clearly shows the transition 
between the linear (Doppler) regime, where 77 3> 1, to the 
quadratic (Dicke) regime, where i/ < 1. The results of 
the numerical integrations throughout the Doppler-Dicke 
transition are well approximated (dashed lines in Fig. [4|) 
by the analytic expression: 

FWHM = 2x (a^^\ , (38) 

where H (x) = e~ x — 1 + x and a 2 = 2/ In 2. The func- 
tion H (x) is usually associated with the velocity self- 
correlation in Brownian motion [l8l |32| , and its extreme 
limits are H(x — > 0) = x 2 /2 and H(x — > 00) = x. 



IV. THE DICKE-DIFFUSION SOLUTION 

In most realistic experiments with EIT, 7 is large 
enough to cause Dicke-narrowing of the two-photon line. 
In principle, this spectrum can be obtained by applying 
the large 7 limit to the general solution of the suscepti- 
bility of the system, Eqs. (|2T|) and (|29[) . However, since 
Eq. ([27|) is somewhat opaque, it seems worthwhile to first 
derive a diffusion-like equation for i?2i(r, t), and then find 
the response of i?3i(r,i). We do this first in the absence 
of fields and then in their presence, analyzing the cases 
of stationary and non-stationary probe. 

A. Diffusion during storage-of-light 

We are interested in describing the dynamics of the 
ground-state populations and coherences, 

< 39 > 

in the absence of electromagnetic fields. This amounts to 
the situation investigated in recent storage-of-light exper- 
iments [30l [3l| , in which a probe beam with a non-trivial 
spatial envelope in the transverse plane was stored in 
an EIT vapor and then retrieved. Following Ref. [33|, 
the storage procedure was described by a linear map- 
ping of the probe field onto the ground-state coherence, 
i?2i (r, t). It was then presumed that the dynamics "in 
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the dark" can be described by the standard diffusion 
process, namely [ifeidiffusion = D\7 2 R G , where D is the 
diffusion coefficient [2f|. Eventually the retrieved probe 
field is obtained from i?2i (r, t) by the reverse linear map- 
ping. 

The dynamics of R21 (r, t) can be described in terms 
of a diffusion-like equation. To this end we return to 
Eqs. (fT3|) . substitute Sli = SI 2 = 0, and apply a procedure 
similar to that of Chapman-Enskog (see e.g. [M])- The 
equation for P21 is 



where 



d ^ . d 

■777 + r 21 + 7 + luJ 2l + V- — 

at or 



P21 = jR 2 i(r,t)F (v) . 



(40) 

Integrating over velocity and using Eqs.© and (fTTj) , we 
have 



^2 + i (^21 - «r 2 i) 

at 



R 2 x(r,t) + ^-3 21 (r,t)=0, (41) 
or 



where 



J ss / (r, t) = J d 3 vvp ss , (r, v, t) 



(42) 



is the current density of the density-matrix. When 7 is 
dominant, we divide Eq. (l40"l) by 7 and get, to zero order 
in I/7, 



D = 



7 



(49) 



is the spatial diffusion coefficient [32|. Substituting this 
result into Eq. (|4"Tj) , we find 

d_ 

at 



R21 (r,t) = DV 2 R 21 (r,t) - (r 21 + iu 21 ) i? 21 (r, i). 



(50) 

Eg. ([50]) describes a spatial diffusion of the coherence, ac- 
companied by a homogenous decay of rate T21 and a 
rotation of rate ^21- A similar derivation can be pre- 
formed for the ground-state populations, i?n and -R22, 
and it results in a similar diffusion equation. This solu- 
tion affirms the theoretical conjectures of Refs. [3fl [3l|. 
In what follows, it is generalized to describe the diffu- 
sion in the presence of the fields, i.e. during slow-light 
propagation. 

It is interesting to note that even in the limit r 2 i — ► 0, 
Eq. (f50|) results in the decay of the total stored- light en- 
ergy. Assuming the stored coherence, R21, is linear in 
the field's amplitude, the total intensity is proportional 
to the integral over |i? 2 i| 2 , which is not conserved and 
always decreases under diffusive spread. A similar obser- 
vation was made in [35j | in the context of stationary light 
pulses that diffuse along the z— direction (the diffusion 
equation therein originates from a different mechanism) . 



p^(r,v,t)=R 21 (r,t)F(v) 



(43) 



B. Diffusion in the presence of fields 



which carries no current. Thus, to find the current we go 
to first order in 1/7, 

p 2 i(r,v,t) = R 21 (r,t)F(v) + -p^(r,^,t), (44) 

7 

substitute it into Eq. (|40p . multiply by Vj and integrate 
over velocity, 



d 3 t 



d d _ 

Vi-^ — h — + r 2 i + 7 + 1U21 



(45) 



1 



, [R 21 (r,t)F(v) + ^{r,v,t) ) = 0. 



Since 



d 3 vv j v i —R 21 (r,t)F{v) = 6 ljV 2 h — R 21 (r,t), (46) 



and 



d 3 vv j p' 2 1 i(r,v,t) = jJ 21lj (r,t), 



(47) 



we find from Eq. ([45| . keeping only leading terms in 7, 



3 21 {r,t) = -D-^R 21 (r,t), 



(48) 



Here we derive the dynamic equations for the envelopes 
of the densities, R 2 i (r, t) and R31 (r, t) , along similar 
lines as above, while considering the interaction with the 
pump and the probe. For brevity, we omit the (r, t) nota- 
tion and denote Sq = q x — q2 . We start from the envelope 
equations, Eqs. lfTT]) . and integrate Eq. (|17b|) over velocity, 



+iSq } J , 1 



^-zA + T 21 ) /7, = /'<>_>( r )./?;; j . 



(51) 

where J ss ' is the envelope of the current densities, defined 
in analogy to Eq. (j4*2"|) . We expand p 2 i as in Eq.(|3I|). 

fax = R21F (v) + (I/7) p ( 2 \\ multiply Eq.(TT7b| by v and 
integrate over velocity. Using Eqs. (|4*6]) and (ITT)) , and 
keeping leading terms in 7, we find 



J21 



dr 



iSq R21 



i J 31 j 

7 



(52) 



where D is defined in Eq. (|49|) . Substituting J 2 i back into 
Eq.([5ID we get 



l-iA + T 21 -D(l + iS q 



R21 



8 \ fWrl- 

/<),(!'!/?,) -;(_+tfqJ._2ilj 31 . 



(53) 
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In order to calculate R31 and J31, we assume in 
Eq. (|17a[) that temporal and spatial changes in the enve- 
lope of the probe are much smaller then the one-photon 
homogenous decoherence rate (Td + 7) and the wave- 
number (qi) , respectively: 



d_ 

dt 



d_ 
dr 



«|6| = |Ai-qrv + z(r d + 7 )|. (54) 



We then formally solve Eq. (|17ap for p 3 i and substitute 

only the dominant part of p 2 \, i.e. p 2 °i = R2iF(v), to 
find 



P3i(r, y,t) 



ijR 31 ~Q 2 (r)R 21 -ni(r,t)no F(v 



(Ai -q r v-H (r d + 7)) 
Integrating Eq. ([55|) over velocity, we get 



(55) 



Rsi = Gi \ijR 3 i - Q 2 (r)R 21 - ni(r,t)noJ , (56) 

where G\ of Eq. (|57|) is the Doppler profile, or 

£31 (r,t) =iK[n 1 (r,t)n + n 2 (r)R 21 {T,t)\ , (57) 

where K = K (Ai) of Eq. ([36|) is the one-photon absorp- 
tion spectrum, which in the Dicke limit of the EIT can 
be considered as a constant near the EIT line. 

Eq. ((53|) and Eq. ([57| form a complete set for R 2 \ and 
R31, only when the term 



0|(r) 



— + is' 

dr ) 



J31 



(58) 



in Eq. (|53|) can be neglected. This term vanishes com- 
pletely in the special case of pump and probe which are 
plane- waves (d/dr =0), collinear and degenerate (<5q = 
0). It can also be neglected whenever |f2 2 | <C 7, as is 
the case in many realistic situations, cf. [21]. Further- 
more, if the latter condition is not satisfied, we can still 
neglect the term (l58l) when both the spatial variations 
(d/dr) and i5q reside in the transverse plane, perpendic- 
ular to qi, since J31 || qi (as can be found by multiply- 
ing Eq. (|55p by v and integrating). Discarding this term 
from Eq. (T5"3")) , and together with Eq. ([57| . we find for i? 21 
a diffusion-like equation, 



-iA + T 21 +K\n 2 {r) 



d 



R21 (r,t) 



(59) 



Eg. (|59)l is the extension of Eg. ([50]) in the presence of 
fields, and it is written in terms of the envelopes. The 
term iSq is responsible for the diffusion across the fields' 
interference pattern, created when qi ^ q 2 . Once we 
solve the diffusion equation for R 2 i(r, t) we substitute it 
in Eg. ([57)) and obtain _R 31 (r, t). We carry this out in the 
next subsection for the case of a plane-wave pump and 
in section [V] for a finite-size pump. 



C. The Dicke-like Absorption Spectrum 

We would like to calculate the susceptibility of the sys- 
tem for the case of a plane- wave pump, f2 2 (r) = £l 2 . 
We Fourier transform in r and t using Eg. (|19p and turn 
Eq. ([5"9]l into 



[i(A + u)-r 21 -K\n 2 \ 2 ]R 21 {k 7 u;) (60) 
= D (Sq + k) 2 R 21 (k, u) + n Kn%nx(k, u), 

and Ea.([5"7|) into 

R 31 (k, w) = iK [Oi (k, w) n + n 2 R 21 (k, w)l . (61) 

Solving Eq.dHDJ) for R 21 (k,u;), substituting in Ea.([6"T]) to 
obtain R 3 i and using Eq. l[2T)]) , we find the susceptibility 
in the diffusion-limit case to be 



X3i(k,u;) = ^An [l-L(k, W )], (62) 
c 



where 

L (k, u) 

and 



i (A + w) - r hom -D(8q + k)' 



horn 



r 21 + K\n 2 \ 



(63) 



(64) 



is the EIT width in the absence of diffusion. 

The absorption of the probe is proportional to Qmx3i ■ 
As explained in section ITTTl in typical cases K is real and 
therefore 



3?mx3i 
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Kn Q [1-KeL (k,w)] 



(65) 



i.e. the well-known EIT absorption spectrum is com- 
posed of the one-photon absorption, gKno/c, and a 
"transparency window" , %ieL (k, u), of Lorentzian shape: 



JieL (k,w) 



K\n 2 \ 2 [r hom + D(Sq + k) 2 } 



(A + o/r + [r 



horn 



D (Sq + k) 



(66) 



Considering L (k, oj) for a given k, as a function of the 
Raman detuning (A + u>), we find the homogenous EIT 
width to be Thorn = ^1 + A^|f2 2 | 2 , broadened by the 
Dicke-EIT width, D (q x q 2 + k) 2 . Eqa.® and ([66]) 
generalize the results of our previous work in Ref . [18J for 
a finite probe in space and time. It includes the power- 
broadening effect, if|f2 2 | 2 that was absent in Ref.[l8|, in 
which the low-contrast approximation was taken. No- 
tice that the term A"|il 2 | 2 replaces the standard power- 
broadening term, |r2 2 | 2 /r c ;, to incorporate the Doppler- 
broadening of the one-photon line. 
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a. 




FIG. 5: The spatial-frequency filter (normalized EIT trans- 
mission), given in Eq. (|66|) , as a function of |kx| = |fex| , for 
different Raman detunings, with lu = 0, <5q = and T2i = 
K\Q 2 \ 2 - The solid-black curve is plotted for A = 0. The 
dashed red and blue curves demonstrate the effect of non- 
zero Raman detuning: a decrease in transparency alongside a 
change in the curvature near k = 0. 



D. Spatial-frequency filter and diffusion-like 
behavior 



When Qmx3i (k, lu) is considered in steady-state (lu — 
0) as a function of k, it acts as a spatial-frequency fil- 
ter for the probe beam. We have in mind a station- 
ary probe beam in the plane z = Z\ with an envelope 
f?i (z\\ x, y), which propagates through the medium to 
the z — z 2 plane. Following Eq. ([M|) . 3mx3i (kj_,cj) de- 
termines the absorption of each spatial- frequency compo- 
nent of Vti (z;k_[_). Since the first term in Eq. (|FJ5)) - the 
one-photon absorption - is constant for all k^ , the filter 
becomes a "transmission" filter with the shape SfteL (kj_) . 
Several examples for this filter are plotted in Fig. [5] with 



2ri2 and <Sq = 0. 



On Raman-resonance, the 
v/r hom /£>, and 



typ 



curve is a Lorentzian of width h 
the maximum overall transmission is achieved. For non- 
zero Raman detuning it obtains a more detailed struc- 
ture - the curvature at k — decreases, becomes zero for 
A = Thom , and turns negative for A > Thorn- 

Figure [6] presents possible outcomes of fii (z2\ x, y), for 
several fli(zi;x 7 y) (representing the absolute value of 
f2i as two-dimensional images). When the features in 
the incident image are large, namely when fii (zi; kj_) 
is confined within k± <C fctyp, the filter is approximately 
quadratic in k± (the central part of the solid- black line in 
Fig. [5]) , which results in pure diffusion in real space. This 
is demonstrated in Figs, [^a) and^b), where the feature- 
size is of the order of 2ir/kt yp . In Fig. [6](c) we demon- 
strate the property of complex diffusion - the phase of 
the left line in fii (zx;x,y) was shifted by n, causing a 
destructive interference between atoms that diffuse to 



(a) 



(b) 
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FIG. 6: Calculated effect of the spatial-frequency filter with 
a plane-wave pump and a finite probe beam: The initial 
pattern, |f2i (zi;x, y) \ (left), and the transmitted pattern, 
|f2i (z2\ x, y) |, after a certain propagation length (right). The 
calculations were done using Eqs. (|30[) . H34|l . (|62|l and (I63|l . 
and the diffraction (k 2 /2/qi) was neglected for clarity. The 
parameters correspond to the black line in Fig. [5](-ftT |il2| 2 = 
Ti2, Sq — 0, A = 0). Images (a)-(c) illustrate regular diffusion 
of real (b) and complex (c) fields. To generate (c), the phase 
of the left line in the incident field was flipped. Images (d) 
and (e) illustrate the effect for smaller features, when larger 
fcx's are pronounce. 



the area between the lines, keeping it dark indefinitely 
[3l| . When smaller features exist and k± extends be- 
yond fctyp a more elaborate behavior occurs. Figure [6{d) 
is a small-scale version of [BJa) and we see, by compar- 
ing Efe) and [f^b) , that the Lorentzian-shaped filter pre- 
serves the sharp edges in the smaller image. This is due 
to the substantial deviation of the Lorentzian filter from 
a quadratic-shaped filter (pure diffusion) for higher fc's. 

A direct measurement of the filter was carried out 
in Refs. [2(J [H|. In these experiments, the trans- 
mission of a stationary (to = 0) on-resonance (A = 0) 
probe beam was measured in two cases: an almost plane- 
wave probe, propagating with a small angular deviation 
from the pump (<5q ^ ; k = 0) , and a divergent probe 
(<5q = ; k^O). In both cases, the wave-number was 
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written as a function of the angular deviation (9), as 
|<5q + k| = q9, and the measured transmission agreed 
with 



IkeL (k) 



Dq 2 6 2 



(67) 



We observe that, in the diverging beam experiment [2l| . 
the outer parts of the beam are absorbed and conse- 
quently the beam radius decreases. This result is counter- 
intuitive from the viewpoint of diffusion, since diffusion 
usually results in spreading, rather than contraction, of 
an initial profile. Nevertheless, in the process of complex 
diffusion, a reduction in size can occur through destruc- 
tive interference - the diffusion wipes out the field in 
regions where the phase pattern exhibits rapid changes. 
The relative phase between a diverging Gaussian probe 
beam and a plane-wave pump beam exhibits rapid vari- 
ations in the transverse direction as the axial distance 
increases, and thus diffusion results in a gradual elimina- 
tion of the outer parts of the beam. 



E. Diffusion of slow-light 

To conclude this section, we analyze the non-stationary 
behavior of slow-light propagation. We consider the case 
of colinear pump and probe, qi || q2 || z, and assume that 
the changes in the probe's envelope along z are much 
smaller than the changes in the transverse plane, such 
that <5q- k <Cfc 2 (or alternatively, take 5q = 0). Denoting 
k = kj_ + k z z 7 with k z <C k±, Eq. ([63|) becomes 



L 



i(A + uj) 



DSq 2 



(68) 



We further take the standard slow-light assumption, that 
the band- width of the probe pulse is fully within the lin- 
ear dispersion regime, i.e. that |A + w| <C \Th om + DSq 2 \. 
More importantly, we assume that Dk 2 ± <C T\ wai + DSq 2 , 
i.e. that the correction to the EIT width resulting from 
the finite size of the probe is much smaller than the EIT 
width of a plane-wave probe. The latter condition pre- 
vents the dispersion of different spatial- frequency compo- 
nents of the envelope, and is the essence of the diffusion 
approximation, allowing us to write L as quadratic in k 



±> 



L 



D5q 2 



1 



i{A + uj)-Dk 2 
Thom + DSq 2 



(69) 



We return to the envelope equations of the probe, 
Eqs. (|3T3)) - (f3Tj) . and use the susceptibility of Eq. 



kz _ 1 + *L _ ii Kno (i _ L ) 



fii(k,a;) = 0. (70) 



1 

c 2qi c 

Substituting L and defining the group-velocity V g as 



gn K 2 \n 2 \ 2 

(r hom + D5q 2 y 



(71) 



we obtain 



ikz - — 



gna 



V„ 



v 9 
1 

c 



K- 



2qi c 
(r h om + D5q 2 - Dk\ 



(72) 

ni(k,£j) = o. 



Returning to the time and space coordinates and assum- 
ing V 9 <c for brevity, we find 



V, 



d d 



9 dz 



dt V2q x 



V a 



D VI + To - iA 



ni(r,t) = 0, 

( ?3 ) 

where Tq = V g gnoK/c — T\ wm — DSq 2 is the on-resonance 
decay rate and is the Laplacian perpendicular to the 
z— axis. Introducing the traveling envelope of the probe 
beam, f^ rav (r,t), as 



fii(r,t) =nf av (r-zV g t,t)e 



(iA-r )t 



(74) 



we find that it undergoes a simple diffusion equation with 
a non-real coefficient: 

J^i rav (r, t) = + D^j Vl&r (r, t) . (75) 

According to Eq. (|75[) . a probe field with an arbitrary 
complex envelope that satisfies the slow-light assump- 
tions will undergo both diffusion, as a result of the atomic 
thermal motion, and optical diffraction. The diffraction 
depends on the actual distance traveled by the beam (due 
to the factor V g /c), while the diffusion depends on the 
time duration. 

An interesting and important example is the propa- 
gation of a paraxial Gaussian beam, such as a Hermite- 
Gaussian or a Laguerre-Gauss mode 36]. These modes 
have the well-known property of being self-similar un- 
der diffraction, i.e. during the propagation through a 
diffractive medium their transverse shape remains un- 
changed up to a length-scale factor. It can be shown from 
Eg. ([75)1 . by utilizing the "complex scaling factor" repre- 
sentation [36j, that Gaussian modes are also self-similar 
under diffusion. For example, the lowest order mode has 
a Gaussian intensity profile, which is known to main- 
tain a Gaussian shape when diffusing. It can further be 
shown from Eq. (|T5[) that the Gaussian modes experience 
a diffusion-induced decay throughout the propagation, as 
discussed at the end of ^IV Al 



V. FINITE PUMP AND PROBE BEAMS 

When both the probe and the pump beams are fi- 
nite, atoms can leave the light beams, evolve "in the 
dark" and diffuse back inside. It was recently demon- 
strated that such a process may result in an EIT line 
much narrower than expected from time-of-flight (TOF) 
broadening and power broadening - a phenomenon de- 
noted as Ramsey narrowing [24| . The line shapes re- 
sulting from TOF broadening and Ramsey narrowing 
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can be described by following the possible atomic paths 
(trajectories), calculating temporal probability functions 
for the atoms' location, and averaging over them (cf. 
[24j). If the time it takes to achieve steady-state with 
the driving field (pumping rate) is comparable to the 
TOF, it is commonly claimed that calculations cannot 
be done with standard steady-state approaches and time- 
dependent solutions have to be used [37J. Nevertheless, 
here we calculate these effects using the steady-state so- 
lution of the diffusion equation in the presence of the 
fields. This is an exact and easier approach that allows 
more elaborate beam geometries to be considered. Note 
that the Ramsey-narrowing experiments are usually done 
with equal pump and probe intensities, while our model is 
for the weak probe regime. However, we expect the main 
attributes of the spectrum to be essentially the same for 
both cases. 

We consider finite probe and pump beams and re- 
strict the discussion to a colinear EIT, qi — q 2 = Sqz. 
We assume that the fields are stationary and overlap in 
their cross sections with a neglected variation along the 
z— direction, 

n 1 (r,*) = niio(r x ) ; n 2 (r) = n 2 w(r±), (76) 

with u>(rj_) the transverse profile of the fields. In the 
diffusion regime we use Eqs. ([5§|) and ([ST)) , which can now 
be written as 

[r - iA + K\n 2 \ 2 \w (rj.) | 2 - DVi] B21 (r±) = 



where 



and 



R 3 i(r±) = iK fi 2 -R2i(r±) + fiino 



w(r ± ) . 



(78) 



where we denoted T = T 2 i + DSq 2 to be the non-power- 
broadened width. In what follows, we solve for R 2 i (rj_) 
and i?3i (rj_) in a specific example and calculate the re- 
sulting absorption spectrum. 



Example: a Stepwise Beam 

We consider a probe and a pump beams with uniform 
intensity and phase within a sheet of thickness 2a in the 
x— direction (one-dimensional stepwise beams): 



w(x,y) 



1 for |x| < a 
for Ixl > a 



(79) 



Eg. ([77)1 can then be written as 



if 



D \ k i-f^) R ^ CM ^ a ) = -noKQ^, 



D kl 



dx 2 



Bai(\x\ >a) = 0, 



(80) 



h = ^(r + K\n 2 \2-iA)/D, 



k 2 = V(r - iA)/D, 



(81) 



and 5Re {fci} > 0. For \x\ < a we expect a solution sym- 
metric in x and for \x\ > a we expect a solution decaying 
for \x\ — > 00. We thus find 



R21 (M < a) = Acosh(kix) 



Kn 
k 2 D 



R21 > a) = 5exp [-fc 2 (|^| - a)] 



(82) 



and the coefficients A and B are obtained from the con- 
tinuity conditions of R 2 \ and -§^R 2 i at \x\ = a: 



Kiln 



kfD cosh(fcia) + {k\/k 2 ) sinh(fcia) 
B = Acosh(fc!a) — 2 fl^i. 



(83) 



For R 3 i(x) we find from Eq.flTS]) that i? 3 i(|x| > a) = 
and 



R3i(\x\ < a) = i£lin K ( 1 



k\D 



(84) 



+ iKAVt 2 cosh (kix) 
The energy absorption at frequency toi is 

P(A) = 2Huj 1 ^- J dzSmj 0^31 (z)} (85) 

and we find 

{„ K 2 \n 2 \ 2 



P(A) = P $lel.K 



T + K\Q 2 \ 2 - iA 



[1 - s D (A); 



where Pq = 2hu>ino and 
tanh(fcia) 



S D (A) 



kia 1 + (ki/k 2 ) taxih(feia) 



(86) 
(87) 



is the correction resulting from the finite size of the beam. 
Figure |7J depicts P(A) and ^eS D (A) for a = lOO/an 
and a —> 00 (plane- wave) with typical parameters. The 
outer part of the finite-beam spectrum (dashed-blue) is 
broadened due to the TOF effect. Ramsey-narrowing is 
apparent in the central part as a cusp-like curve. The 
cusp is narrower than the power-broadened Lorentzian, 
but it is nevertheless limited by the width T = T 2 i + 
DSq 2 . 

A similar calculation can be done for the two- 
dimensional analogue of the stepwise sheet: a stepwise 
cylindrical profile with w (r < a) = 1 and w (r > a) = 0, 
where r 2 — x 2 + y 2 . Solving Eqs. (|77|) and (|75|) in cylin- 
drical symmetry, one finds that the energy absorption 
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FIG. 7: Normalized transmission for a plane-wave and a 
finite-sized beam (ID and 2D), demonstrating TOF broad- 
ening and Ramsey-narrowing. The inset depicts the correc- 
tion to the spectrum, resulting from the finiteness of beams, 
Ke [1 - S D (A)], for the ID case [Eq.fg?])]. The parameters 
are: T = 100 Hz, K\n 2 \ 2 = 2 KHz, D = 10 cm 2 /sec and 
a — 100 fim. The choice of A'lf^l 2 3> T makes the narrowing 
effect more obvious. 



spectrum, -P(A), has the exact same form of Eq. (|86|) . 
with the correction term being 

• (88) 



I (ha) ki K (k 2 a) 



h (kia) k 2 K\ (/c 2 a) 



S D (A) = — 

Here, k\p are as defined in Eqs.JH]) and Iq (x) and Kq (x) 
are modified Bessel functions. A comparison between 
the ID and the 2D spectra is given in Fig. [7] In the 
2D case the TOF effect is substantial while the Ramsey- 
narrowing is reduced. The latter can be attributed to the 
fact that on average less atoms return to the beam in the 
2D geometry. 



VI. CONCLUSIONS 

We have presented a model for EIT that incorporates 
thermal atomic motion by introducing the density-matrix 
distribution in space and velocity along with a Boltzmann 
relaxation term. The model describes a range of motional 
phenomena, such as Dicke narrowing, Ramsey narrowing 
and diffusion during storage of light, which have been 
analyzed in the past in different independent studies. In 
the absence of electromagnetic fields, the model reduces 
to pure diffusion of the ground-state coherence and pop- 
ulation, in agreement with recent storage of light experi- 
ments. In the presence of a pump and a weak probe fields, 
the solution for the spatial distribution of the atomic co- 
herence and the fields can be obtained. 



For the case of a plane-wave pump, the temporal and 
spatial frequency components of the incoming probe beam 
form the natural basis for the problem. The general 
solution, in terms of the complex susceptibility of the 
medium, X3i(k, lu), exhibits the Doppler-Dicke transi- 
tion for both the one-photon and the two-photon absorp- 
tion spectra. For the realistic regime, in the presence 
of a buffer-gas, when the one-photon line is Doppler- 
broadened and the two-photon line is Dicke-narrowed, 
an explicit expression for the EIT transmission, L(k, lu), 
is derived. £(k, lu) yields the EIT absorption spectrum, 
for any given k, including the power-broadening effect 
and the Dicke width. Moreover, for any given lu, L(k, lu) 
serves as a spatial-frequency filter, generally diminishing 
high k values. We explain this diminution by the diffu- 
sion of atoms across the pump-probe interference pattern, 
of wave- vector k + qi — q 2 . On Raman-resonance, and 
when the spatial features of the incoming probe beam 
are large, the probe's envelope undergoes a diffusion-like 
dynamics. For smaller features or for non-zero Raman- 
detuning, a more elaborate behavior takes place. We note 
that a unique result is obtained for probe beams with a 
single value of |kj_|, usually referred to as non-diffracting, 
e.g. Bessel beams [38j . These will not be distorted by 
L(k, lu) and hence will not spread due to neither diffusion 
nor diffraction. 

Since the complex amplitude of the probe's envelope 
diffuses, interference phenomena occur. For example, 
destructive interference between adjacent features that 
are opposite in phase, maintains the dark area between 
them. This also explains why adjacent rings in the 
Bessel beam remain separated indefinitely. Furthermore, 
when the optical diffraction is taken into account, the 
effective diffusion coefficient becomes a complex number 
(D + iVg/2/qi), with the group-velocity determining the 
ratio between the real (pure-diffusion) and the imaginary 
(diffraction) parts. Such complex diffusion can possibly 
be useful for all-optical image processing, such as image 
enhancement, denoising and edge-detection [39| . 

Ramsey narrowing occurs when the pump's cross- 
section is finite, and atoms that re-enter the beam from 
outside, less affected by power-broadening, contribute to 
the spectrum. Our model gives simple and analytic re- 
sults for the Ramsey-narrowed spectrum by solving a 
diffusion equation with spatially dependent decay coef- 
ficients and sources. It is somewhat surprising that our 
steady-state approach is able to capture this effect, for 
any pumping and transit rates, so that one is not required 
to average over atomic trajectories. Utilizing our model, 
the spectrum for any detailed geometry can readily be 
obtained. 

The theory presented here may contribute to the anal- 
ysis of decoherence in collective light memories, for which 
it was shown that the decay rate is proportional to the 
single-atom diffusion rate [4fj|. The model can poten- 
tially be extended to include the main ground-state deco- 
herence mechanisms in vapor: non-coherence-preserving 
collisions, namely spin-exchange collisions, and wall col- 
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lisions. Such extensions may aid in developing vapor EIT 
schemes with narrower lines. 
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